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Abstract 

Two permutations of [n] are comparable in the Bruhat order if one is closer, in a natural way, to the 
identity permutation, 1 2 • ■ ■ n, than the other. We show that the number of comparable pairs is of order 
(n!) 2 /n 2 at most, and (n!) 2 (0.708)™ at least. For the weak Bruhat order, the corresponding bounds are 
(n!) 2 (0.362) n and (n!) 2 FJiLi (H (*) A)> where H (i) := X^-=i Vi- ^ n light of numerical experiments, we 
conjecture that for each order the upper bound is qualitatively close to the actual number of comparable 
pairs. 



Introduction 

Let n > 1 be an integer. Two permutations of [n] := {1, . . . ,n} are comparable in the Bruhat order if 
one can be obtained from the other by a sequence of transpositions of pairs of elements forming an inversion. 
Here is a precise definition of the (strong) Bruhat order on the set of permutations S n (see Stanley p. 
172, ex. 75. a.], Humphreys |12l p. 119]). If lu = lu (1) • • • to (n) € S n , then a reduction of lu is a permutation 
obtained from to by interchanging some tu (i) with some lu (j) provided i < j and lu (i) > lu (j). We say that 
7r < er in the Bruhat order if there is a chain a = u>\ — *■ u>2 — * • • • — ► u> s = it, where each io t is a reduction of 
LOt—x- The number of inversions in Lo t strictly decreases with t. Indeed, one can show that if u>2 is a reduction 
of uji via the interchange toi(i) <-> Wi(i), i < j, then 

inv(wi) = inv(cj2) + 27V(wi) + 1, 
N(wi) :=\{k : i<k< j, > oj^k) > wi(j)}|; 

here inv(wi), say, is the number of inversions in lo\ (see Bjorner and Brcnti [2]). Figure 1 below illustrates 
this poset on S3 and S4. 

The definition of the Bruhat order is very transparent, and yet deciding for given 7r, a whether 7r < a 
is computationally difficult, even for smallish n. Fortunately, the historically first definition of Bruhat 
comparability, due to Ehrcsmann (1934) 9 , is actually an efficient algorithm to check if ix < a (see also 
Knuth |T2], El)- The Ehresmann "tableau criterion" states that 7r < a if and only if iTi.j < ai.j for all 
l<i<j"<jT, — 1, where TTij and o~i .j are the i-th entry in the increasing rearrangement of it (1) , . . . , it (j) 
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and of a (1) , . . . , er (j). These arrangements form two staircase tableaux, hence the term "tableau criterion" . 
For example, 41523 > 21534 is verified by element-wise comparisons of the two tableaux 

1 2 4 5 1 2 3 5 

14 5 12 5 

14 12 
4 2 

Then, Deodhar |7] extended the Bruhat order notion to other Coxeter groups. It is well-known that Ehres- 
mann's criterion is equivalent to the (0, l)-matrix criterion. It involves comparing the number of l's contained 
in certain submatrices of the (0, l)-permutation matrices representing 7r and a (see Bona [5], Later, 
Bjorner and Brcnti |3] were able to improve on the result of [5], giving a tableau criterion that requires 
fewer operations. Very recently, Drake, Gerrish and Skandera |5| have found two new comparability crite- 
ria, involving totally nonnegative polynomials and the Schur functions respectively. We are aware of other 
criteria (see Bjorner 0], Fulton [111 pp. 173-177], Lascoux and Schutzenberger 0), but we found the 
(0, l)-matrix and Ehresmann criteria most amenable to probabilistic study. 
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Figure 1: The Bruhat order on S3 and S4,. 



The (0, l)-matrix criterion for Bruhat order on S n says that for 7r, a € S n , tt < a if and only if for all 
h j < n i the number of 7T (1) , . . . , 7T (i) that are at most j exceeds (or equals) the number of a (1) , . . . , a (i) 
that are at most j (see |S] for this version). It is referred to as the (0, l)-matrix criterion because of the 
following recasting of this condition: let M (tt), M(o~) be the permutation matrices corresponding to tt, a, 
so that for instance the (i, j)-entry of M(ir) is 1 if tt (j) = i and otherwise. Here, we are labeling columns 
1, 2, ... ,71 when reading from left to right, and rows are labeled 1, 2, ... ,71 when reading from bottom to 
top so that this interpretation is like placing ones at points (i, tt (i)) of the n x n integer lattice and zeroes 
elsewhere. Denoting submatrices of M (•) corresponding to rows / and columns J by M (-) 7 J; this criterion 
says that tt < a if and only if for all i, j < n, the number of ones in M Mm ui is at least the number of ones 
in M(cr)j i j (see jS] for this version). 

An effective way of visualizing this criterion is to imagine the matrices M (tt) and M (a) as being superim- 
posed on one another into a single matrix, M (it, a), with the ones for M (it) represented by x 's ("crosses"), 
the ones for M (a) by o's ("balls") and the zeroes for both by empty entries. Note that some entries of 
M (tt, a) may be occupied by both a cross and a ball. Then the (0, l)-matrix criterion says that 7r < a if 
and only if every southwest submatrix of M (tt, a) contains at least as many crosses as balls. Here, in the 
notation above, a southwest submatrix is a submatrix M (ir, c)^ ^ of M (tt, a) for some i, j < n. It is clear 
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that we could also check tt < a by checking that crosses are at least as numerous as balls in every northeast 
submatrix of M(ir,a). Likewise, tt < a if and only if balls are at least as numerous as crosses in every 
northwest submatrix of M (tt, a), or similarly balls are at least as numerous as crosses in every southeast 
submatrix of M (tt, a). Parts of all four of these equivalent conditions will be used in our proofs. As a quick 
example, with tt = 21534 and a — 41523, tt < a is checked by examining southwest submatrices of M (tt, a) 
in Figure 2 below. Also, the superimposing of M (tt) with M (a) to form M (tt, a) is illustrated in this figure. 
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Figure 2: Superimposing M(tt) and M(a) to form M(ir,a). 



In this note, we use the (0, l)-matrix and the Ehresmann criteria to obtain upper and lower bounds for the 
number of pairs (tt, a) with tt < a. 



Theorem 1 Let n > 1 be an integer, and let tt, a <E S n be selected independently and uniformly at random. 
Then there exist universal constants c\ , c 2 > such that 

ci (0.708)™ < P (tt < a) < c 2 /n 2 . 



Equivalently, the number of pairs (7r,cr) with tt < a is sandwiched between ci(0.708)"(n!) 2 and c 2 n~ 2 (nl) 2 . 
The lower bound follows from a sufficient condition derived from the (0, l)-matrix criterion, and a computer- 
aided tabulation of an attendant function of a smallish integer argument. Empirical estimates based on 
generating pairs of random permutations suggest that P (tt < a) is of order n~( 2+s \ for S close to 0.5 from 
above. So apparently it is the upper bound which comes close to the true proportion P (tt < a). It is certain 
that the constant 0.708 can be further improved, but we do not know if our method could be extended to 
deliver a lower bound (1 — o(l))™. A lower bound n~ a , a qualitative match of the upper bound, seems out 
of sight presently. 

Then we turn to the modified order on S n , the weak Bruhat order Here tt ^ a if there is a chain 

a = u>i — > u>2 — > • • • — * u s = TTj where each u> t is a simple reduction of u>t-iy i-e. obtained from u>t-i by 
transposing two adjacent elements u)t-i(i), u>t-i(i + 1) with u)t-i(i) > cJt-i(i + !)• Since at each step the 
number of inversions decreases by 1, all chains connecting a and tt have the same length. 



Theorem 2 Let P* := P (tt <cr). Then P* is submultiplicative, i.e. P* 1+n2 < P^P^- Consequently there 
exists p = lim y/P*. Furthermore, there exists an absolute constant c > such that 
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f[(H(i)/i)<P:<c(0.362) n , 

i=l 

where H (i) := X)}=i Vi- Consequently, p < 0.362. 

The proof of the upper bound is parallel to that of Theorem ^ lower bound, while the lower bound follows 
from an inversion set criterion inspired by discussion of the weak Bruhat order by Berge P^. Empirical 
estimates indicate that p is close to 0.3. So here too the upper bound seems to be qualitatively close to 
the actual probability P*. And our lower bound, though superior to the trivial bound 1/n!, is decreasing 
superexponentially fast with n, which makes us believe that there ought to be a way to vastly improve it. 
Paradoxically, it is the lower bound that required a deeper combinatorial insight. We show first that the 
number of 7r's below (or equal to) a equals e('P), the total number of linear extensions of V = Tier), the 
poset induced by a. (The important notion of Via) was brought to our attention by Sergey Fomin |10|.1 
And then we prove that for each poset V of cardinality n, 

e(V)>n\/l[d(i), 

where d (i) := | {j G V : j < i in V} |. The final step is based on the independence of sequential ranks in the 
uniform permutation. Reducing the gap between the bounds will probably require a better understanding 
of V{o) for a typical permutation a. 

In conclusion we mention two papers, |18| and |19| . where the "probability-of-comparability" problems 
were solved for the poset of integer partitions of n under dominance order, and for the poset of set partitions 
of [n] ordered by refinement. 

Proof of Theorem [T], upper bound. 

We need to show that 

P(n < a) = 0{n- 2 ) . 

The argument is based on the (0, l)-matrix criterion. We assume that n is even. Only minor modifications 
are necessary for n odd. 

STEP 1. The (0, l)-matrix criterion requires that a set of n 2 conditions are met. The challenge is to select 
a subset of those conditions which meets two conflicting demands. It has to be sufficiently simple so that 
we can compute (estimate) the probability that the random pair (tt, a) satisfies all the chosen conditions. 
On the other hand, collectively these conditions need to be quite stringent for this probability to be o(l). In 
our first advance we were able (via Ehresmann's criterion) to get a bound 0(n -1 / 2 ) by using about 2n x l 2 
conditions. We are about to describe a set of 2n conditions that does the job. 

Let us split the matrices M (it, a), M (it) and M (a) into 4 submatrices of equal size n/2 x n/2 - the 
southwest, northeast, northwest and southeast corners, denoting them M sw (•), M ne (•), M nw (•) and M se (•) 
respectively. In the southwest corner M sw (tt,o~), we restrict our attention to southwest submatrices of the 
form i x n/2, i = 1, . . . , n/2. If n < er, then as we read off rows of M sw (-7T, a) from bottom to top keeping 
track of the total number of balls and crosses encountered thus far, at any intermediate point we must have 
at least as many crosses as balls. Let us denote the set of pairs (ir, a) such that this occurs by £ sw . We draw 
analogous conclusions for the northeast corner, reading rows from top to bottom, and we denote by £ ne the 
set of pairs (it, a) satisfying this condition. 

Similarly, we can read columns from left to right in the northwest corner, and here we must always have 
at least as many balls as crosses. Denote the set of these pairs (ir, a) by £ nw . The same condition holds for 
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the southeast corner when we read columns from right to left. Denote the set of these pairs (n, a) by £ se . 
Letting £ denote the set of pairs (it, cr) satisfying all four of the conditions above, we get 

< c} c £ = £ sw n £ ne n £ nw n£ se . 

Pairs of permutations in £ satisfy 2n of the n 2 conditions required by the (0, l)-matrix criterion. And unlike 
the set {it < a}, we are able to compute \£\, and to show that P{£) = (n\)~ 2 \£ \ = 0(n~ 2 ). Figure 3 below 
is a graphical visualization of the reading-off process that generates the restrictions defining the set £. 



M(n,o) 




Figure 3: Finding a necessary condition for n < a. 



If a row (column) of a submatrix M(ir)i t j (M(a)ij resp.) contains a marked entry, we say that it 
supports the submatrix. Clearly the number of supporting rows (columns) equals the number of marked 
entries in M(n)i t j (M(a)i.j resp.). Now, given n, cr, let Mi = Mi(tt), M 2 = M 2 (a) denote the total 
number of rows that support M sw {tt) and M sw (o) respectively. Then M nw (ir), M nw {o~) are supported by 
M 3 = n/2 — Mi columns and by M 4 = n/2 — M 2 columns respectively. The same holds for the southeastern 
corners of M(tt) and M(a). Obviously the northeastern submatrices of M(ir) and M(a) are supported by 
Mi rows and Mi rows respectively. Then we have 

P(£)= £ P(£nA( mi ,m 2 )), (1) 

mi ,m2 

i(mi,m 2 ) := {(tt, a) : Mi = m u M 2 = m 2 } . 
Clearly £ n A (mi, m 2 ) = if mi < m 2 . We claim that, for mi > m 2 , 

4 nlAt) 



P(£ , n^(mi,m 2 )) 



(mi -m 2 + l)(n/2+l) 



(m 3 + l)(mi + 1) 



\n/2) 



Here and below m 3 := n/2 — mi and m 4 := n/2 — m 2 stand for generic values of M 3 and M 4 in the event 
A(m 1 ,m 2 ). 

To prove (f), let us count the number of pairs (tt, a) in £ n ^4(mi,m2). First consider the southwest 
corner, M sw (ir, a). Introduce L\ = L\ (w, a), the number of rows supporting both M sw (ir) and M sw (<j). So 
L\ is the number of rows in the southwest corner M sw (tt, a) containing both a cross and a ball. Suppose 
that we are on the event {Li — We choose l\ rows to support both M sw (tt) and M sw (a) from the n/2 
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first rows. Then, we choose (mi — I \ + 777,2 — £1) more rows from the remaining (77/2 — £\) rows. Each of 
these secondary rows is to support either M sw (ir) or M sw (o~), but not both. This step can be done in 

fn/2\ f n/2 - £ x \ 
\h )\m 1 -£ 1 + m 2 -£j 

ways. Next, we partition the set of (mi — £ x + m 2 — l\) secondary rows into two row subsets of cardinality 
(mi ~ £%) an( i ( m 2~^i) that will support M sw (tt) and M sw (a), accompanying the £ x primary rows supporting 
both submatrices. We can visualize each of the resulting row selections as a subsequence of (1, ... , ?7/2) which 
is a disjoint union of two subsequences, with £\ and (mi — £\ + m 2 — £\) elements respectively, with each 
element of the subsequence marked by a cross and/or ball dependent on which submatrix is supported by 
this element, i.e. the row represented by this element. The condition £ sw is equivalent to the restriction: 
moving along the subsequence from left to right, at each point the number of crosses is not to fall below the 
number of balls. Obviously, no double-marked element can cause violation of this condition. Thus, our task 
is reduced to determination of the number of (mi — £\ + m 2 — £\)-\ortg sequences of mi — £\ crosses and 
m 2 — £1 balls such that at no point the number of crosses is strictly less than the number of balls. By the 
classic ballot theorem (see Takacs 22 pp. 2-7]), the total number of such sequences equals 

(mi — l\ + 1) — (m 2 — £1) /mi — £1 + ni 2 — t\ + 1\ _ mi — m 2 + 1 (m\ — £\ + ni2 — £\ 
(mi - Ii + 1) + (m 2 - £1) V mi- £i + l / mi-^i + 1 \ mi - £1 

The second binomial coefficient is the total number of (mi — £1 + 7772 — ^i)-long sequences of (mi — £1) crosses 
and (7772 — £1) balls. So the second fraction is the probability that the sequence chosen uniformly at random 
among all such sequences meets the ballot theorem condition. The total number of ways to designate the 
rows supporting M sw (jr) and M sw (o~), subject to the condition £ sw , is the product of two counts, namely 



fn/2\/ n/2-£i \ (m x - £ x + m 2 - £A mi - m 2 + 1 
V £1 J \mi -£i+m 2 -£i)\ mi - 4 / mi - £ x + 1 

_ mi — m 2 + 1 / n/2\ ( m 2 \ ( n/2 — m 2 + 1\ 

~ n/2 -m 2 + i U2/ \4y V m i ~ £ i + 1 / 

Summing this last expression over all £\ < 7772 , we obtain 



mi - m 2 + 1 (nj2\ ^ /m 2 \ /n/2 - m 2 + \\ _ m 1 - m 2 + 1 fn/2\ (n/2 + 1\ 
n/2 - m 2 + 1 \m 2 J Z--' \£\ ) V mi - £ x + 1 / n/2 - m 2 + 1 V m 2 / Vmi + 1 J 

x y £i<ni2 x/x // \/\ / 

_ (mi - m 2 + l)(n/2 + 1) /7i/2\ /t7/2\ 
(n/2 — m 2 + l)(mi + 1) \mi J \m 2 J 

Here, in the first equality, we have used the binomial theorem. The product of the two binomial coefficients 
in the final count is the total number of row selections from the first n/2 rows, mi to contain crosses and 
7772 to contain balls. So the fraction preceding these two binomial factors is the probability that a particular 
row selection chosen uniformly at random from all such row selections satisfies our ballot condition "crosses 
never fall below balls" . Equivalently, by the very derivation, the expression J2J is the total number of paths 
(X (t) ,Y (t)) < t <n/2 on the square lattice connecting (0,0) and (mi,m 2 ) such that X(t + 1) — X(t), Y(t + 
1) - Y(t) G {0, 1}, and X{t) > Y(t) for every t. (To be sure, if X(t + 1) - X(t) = 1 and Y(t + 1) - Y{t) = 1, 
the corresponding move is a combination of horizontal and vertical unit moves.) 

Likewise, we consider the northeast corner, M„ e (7r, a). We introduce L 2 = L 2 (ir, a), the number of rows 
in M ne (ir,a) containing both a cross and a ball. By initially restricting to the event {L 2 — £ 2 }, then later 
summing over all £ 2 < m 2 , we obtain another factor Analogously, a third and fourth factor J3J| comes 
from considering columns in the northwest and southeast corners, M nw (n,a) and M se (jr,a). Importantly, 
the row selections for the southwest and the northeast submatrices do not interfere with the column selections 
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for the northwest and the southeast corners. So by multiplying these four factors (0 we obtain the total 
number of row and column selections on the event A(mi, 7712) subject to all four restrictions defining £\ 

Once such a row-column selection has been made, we have determined which rows and columns support 
the four submatrices of M(tt) and M(a). Consider, for instance, the southwest corner of M(ir). We have 
selected mi rows (from the first n/2 rows) supporting M sw (tt), and we have selected 7713 columns (from 
the first n/2 columns) supporting M nw {%). Then it is the remaining n/2 — 777,3 ~ m i columns that support 
M sw (tt) . The number of ways to match these m\ rows and m\ columns, thus to determine M sw (it) completely, 
is 7774!. The northeast corner contributes another 7771!, while each of the two other corners contributes 7773!, 
whence the overall matching factor is (milms!) 2 . The matching factor for a is (m,2^m,4\) 2 . Multiplying the 
number of admissible row-column selections by the resulting Ili=i( m «') 2 anc ^ dividing by (n!) 2 , we obtain 



P(£nA(m 1 ,m 2 )) 



(mi - m 2 + 1) (77/2 + 1) /n/2 
l)Ul 



77/2 
m 2 



nti(^o 2 

(n!) 2 



(n/2 — 7772 + l)(777l 

which is equivalent to (f). Figure 4 below is a graphical explanation of this matching factor. In it, we 
show the matrix M(tt) in a case when in the southwest and the northeast squares tt is supported by the 
bottom m(= mi) and the top m rows respectively; likewise, in the northwest and the southeast squares tt is 
supported by the n/2 — m leftmost and the n/2 — m rightmost columns respectively 
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n/2 - m 

Figure 4: Selection of first m = mi (n/2 — m resp.) rows (columns resp.) in corners to support M(n) 
STEP 2. Let us show that and (f) imply 



P{£)<E 
First, Mi and M2 are independent with 



(Mi - M 2 + l) 4 (n/2 + l) 4 
(77/2 - M 2 + l) 4 (Mi + l) 4 



(±) 



P(M i =rrH) = f¥ v , i = l,2. 
{n/2) 
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Indeed, Mi obviously equals the cardinality of the intersection with [n/2] of a uniformly random subset of 
size n/2 from [n], which directly implies these formulas. Thus, each Mi has the hypergeometric distribution 
with parameters n/2, n/2, n/2; in other words, Mi has the same distribution as the number of red balls in 
a uniformly random sample of n/2 balls from an urn containing n/2 red balls and n/2 white balls. By the 
independence of M± and M 2 , we obtain 



P(Mi = m 1 ,M 2 = m 2 ) 
It remains to observe that Q and (f ) imply 



/n/2\ 2 rn/2y 

(;; 2 ) 2 



r>/ c\ (mi-m 2 + l) 4 (n/2+l) 4 

= 2^ TIT, 7~^A~, ' P ( Ml = mi > M 2 = TO2) 



m 1 > rn o 



{n/2 -m 2 + if (mi + 1)" 



< 2^ TTTT7 ^4 • P(Mi = mi,M 2 = m 2 ) 



( n / 2 - ™2 + ir k + 1) 



= e 



(Mi-M 2 + l) 4 (n/2 + l) 4 



(n/2 - M 2 + l) 4 (Mi + l) 4 



and (|) is proved. 



STEP 3. The advantage of (J) is that it allows us to use probabilistic tools exploiting the independence 
of the random variables M% and M 2 . Typically the Mj's are close to n/4, while |Mi — M 2 | is of order n 1 / 2 
at most. So, in view of (J) we expect that P(£) = 0(n~ 2 ). 

We now make this argument rigorous. First of all, by the "sample- from-urn" interpretation of Mj, 



I n-l \ 
„,„, n \ n /2-l 1 . 



2 (»%) 



(3) 



Then (see Janson et al. ^31) the probability generating function of Mi is dominated by that of Bin(n, 1/4), 
and consequently for each t > we have 

P {\Mi - n/4| >t) = (exp (-4i 2 /n)) . 
Hence, setting t — n 2 / 3 we see that 



P (n/4 - n 2/3 < Mi < n/4 + n 2/3 ) > 1 
for some absolute constant c > 0. Introduce the event 

2 



,1/3 



H{n/4- 



- n 2/3 <M„ <n 



1=1 



/4 + n 2 / 3 } 



Combining the estimates for Mi, we see that for some absolute constant ci > 0, 

P(A n ) > l-e- Cl ™ 1/3 . 

Now the random variable in (|), call it X n , is bounded by 1, and on the event A n , within a factor of 
l + Oin- 1 / 3 ), 



X„ = I ^ ) (Mi - M 2 + l) 4 (n/2 + l) 4 
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Therefore 



P {£) < (~) (n/2 + l) 4 E [{Mi - Mi + l) 4 ] + ( e - Cl " 1/3 ) . 
It remains to prove that this expected value is 0(n 2 ). Introduce Mi = Mi — £7[Mj], i — 1,2. Then 

(Mi - M 2 + l) 4 = (Mi - M 2 + I) 4 < 27(M? + m\ + 1), 

as 

{a + b + cf < 3(a 2 + 6 2 + c 2 ). 

We now demonstrate that E^M^] — 0(n 2 ). To reduce computations to the minimum we use the fact that, 
as a special instance of the hypergeometrically distributed random variable, Mi has the same distribution as 
the sum of n/2 independent Bernoulli variables Yj S {0, 1} (see Vatutin and Mikhailov alternatively ^3 
p. 30]). So denoting pj = P(Yj = 1), qj = P(Yj = 0), we get a product formula for the moment generating 
function of Mi\ 



E 



J/, 



n/2 n/2 

3=1 3=1 

"/ 2 / 2 3 4 

II ( 1 + T p i q i + 3fft'?*fe ~ ft') + 7r^' 9 J 4 + p 3%) + °^ 5 ) 



3=1 



9 n/2 o n/2 

1 + tLw* + ¥L Pi* ~ p 3 > 



3=1 



3 = 1 



/2 



31 #32 



3 = 1 



+ o(« 5 ). 



In particular, 



and consequently 



Var[Mi] = E MA =Y^ 



i/2 



3 = 1 



n/2 

E [Mi] =6 X {Ph<lh){Ph<lh) + J2(Pj1j + PjQj) 

3l#32 3 = 1 

< 6Var 2 [M,] + Var[M,]. 
It is left to show that Var[M;] = O(n). Extending the computation in ©, 



o 

71 

= E 2 [M t ] + 0{n). 
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Therefore 

Var[Af,] = E[M l {M l - 1)] + E[M t ] - E 2 [Mi] = 0(n). 
This completes the proof of Theorem ^ (upper bound). D 

Proof of Theorem [T], lower bound. 

We will show that for each e > 0, 

P(ir<a)=Sl((a-e) n ), 

where 



n/25497938851324213335 „ „„„„„ 

a = \ - = 0.70879 .... 

V 22! 

Introduce 7r* (er* resp.), the permutation 7r [a resp.) with the element n deleted. More generally, for k < n, 
7T (<r k * resp.) is the permutation of [n — k] obtained from ir (a resp.) by deletion of the k largest elements, 
n, n — 1, . . . , n — k + 1. The key to the proof is the following 

Lemma 1 Let k 6 [n]. If every northeastern submatrix of M(tt,o~) with at most k rows contains at least as 
many crosses as balls, and ir k * < a k * , then ir < a . 

Before proceeding with the proof, we introduce one more bit of notation. Let G n be the empty n x n 
grid depicted in the M(-)'s of Figure 2 above. Figure 5 below is a depiction of G5 and an emboldened 
northeastern-corner 3x4 subgrid of it, denoted by C. 



Figure 5: G5 and an emboldened subgrid C. 

If C is any subgrid of G n , then M(- \ C) denotes the submatrix of M(-) that "sits" on C. To repeat, the 
(0, 1) -matrix criterion says that n < a if and only if for each northeastern-corner subgrid C of G n , we have 
at least as many crosses as balls in M(%, cr\C). 

Proof. By the assumption, it suffices to show that the balls do not outnumber crosses in M(tt, cr I C) for 
every such subgrid C with strictly more than k rows. Consider any such C . Let denote the subgrid 
formed by the top k rows of C. Given a submatrix A of M(ir) (of M{a) resp.), let \A\ denote the number 
of columns in A with a cross (a ball resp.). We need to show 

|M(7r|C)|>|M(a|C)|. 

By the assumption, we have \M(tt\C^)\ > \M(a\C^)\. Write \M(n \ C^)\ = \M{a | C<»)| + A, A > 0. We 
now delete the top k rows from M(tt), M(o~) together with the k columns that contain the top k crosses in the 
case of M(tt) and the k columns that contain the top k balls in the case of M(a). This produces the matrices 
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M(7r fc *) and M(a k *). In either case, we obtain the grid G„_fc together with a new northeastern subgrid: 
C(w k *) in the case of M(ir) and C(a k *) in the case of M(a). Figure 6 below is a graphical visualization 
of this deletion process in the special case tt — 12534, a — 45132, k — 2 and C the 3x4 northeastern 
subgrid of G 5 . We have emboldened C in M(tt), M(a), and the resulting C(tt 2 *), C{a 2 *) in M (tt 2 *), M (a 2 *) 
respectively. 
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Figure 6: Deletion of 2 largest elements of tt, <j, and its affect on C. 



Since we delete more columns in the case of tt than cr, note that C(7r fe *) C C(er fe *) as northeastern subgrids 
of G„_fe. In fact, these grids have the same number of rows, but C(7T fc *) has A fewer columns. Hence, as 
TT k * < a k * , we have 

\M{ir k * | C{ir k *))\ > \M{a k * | C{ii k *))\ > \M(a k * \ C(a k *))\ - A. 

So 

\M{tt \C)\ = \M(tt I C (fc) )| + \M{TT k * | C(ir k *))\ 

= \M(a | G (fc) )| + A + \M(ir k * \ C(n k *))\ 
> \M(a | G (fe) )| + \M(a k * | C(a k *))\ 
= \M(ct\C)\, 

which proves the lemma. D 

For each k < n, let £ n ,k denote the event "every northeast submatrix of the top k rows has at least as many 
crosses as balls". Then by Lemma ^ 

{tt < a) D S n , k n {n k * <a k *}. 
Now the events £ n ^ and {7r fe * < a k *} are independent! So we get 

P(7r<a) >P{S n , k )P{TT k * <a k *). (4) 

For the permutation tt (a resp.) introduce £i(Tr) = 7r _1 (i) (^i(c) = cr _1 (i) resp.), the index of a column 
that contains a cross (a ball resp.) at the intersection with row i. In terms of the £i(-)'s, £ n: k is the event: 
for each integer j < k and m < n, the number of £ n (Tr), £„_i(7r), . . . , £ n ^j + i{iT) that are m at least is more 
than or equal to the number of £ n {cr) 7 £ re _i(cr), . . . ,£ n -j + i(a) that are m at least. We could have replaced 
an integer m < n with a real number, which means that 
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£ n , k = {{w,a) : {£{*),£(<?)) eC k }, 

for some cone-shaped (Borel) set C k C K 2fc ; here £(ir) = {£ n - i+ i(n)}i<i< k) £{&) = {i n -i+i(v)}i<i<k- 

Our task is to estimate sharply P (£ n ,k) for a fixed k, and n — > oo. Observe first that £(77) and £(cr) are 
independent, and each uniformly distributed. For instance 

f(4W = il,--- A-fc+iM =ife) = tt") 1 <ii ^ ••• ^ jfc < n, 

where (n)/s = n(n— 1) ■ ■ • (n— fc+1). Since (n)fc ~ n k as n — > 00, ■ ■ ■ A-HiW are almost independent 

[n]-uniforms for large n, and fixed fc. Let us make this asymptotic reduction rigorous. Let {/ be a uniform- 
[0, 1] random variable, and let Ui,...,U n be independent copies of U. Then each \nUi\ is uniform on [n], 
and it is easy to show that 

P (\nU{\ = i u . . . , \nUk\ = i k I \nU{\ + ■ ■ ■ ? \nU k ]) = 

(n) k 

In other words, {£n-i+i(^)}i<i<k has the same distribution as [nU] := { \nUi\ }\<i<k conditioned on 
the event A n , k — {\nU{\ 7^ ■ ■ • 7^ |~n[/fe"|}. Analogously {£ n -i+i{o')}i<i<k is distributed as fnV] := 
{|~nV;]}i<i<fc conditioned on S n) fe = {[nT^i] 7^ ■ • • 7^ [nVfc]}, where V\,...,V k are independent [0,1]- 
uniforms, independent of Z7i, . . . , We will need yet another event T> n ^ on which 



minjmin \Ui — Uj\, min | — Vj\, min \Ui — Vj\} > l/n. 

Clearly on V nyk 



(\nU],\nV])eC k <^(U,V)eC k ; 
here U := {Ui}i<i< k , V := {Vi}i<i<fc. In addition V n _ k C „4 nfc n B ni fe, and 

P(X£ |fc ) < 2k 2 P(\U 1 - £7 2 | < l/n) < 4k 2 /n. 

Therefore 



P(f n , fe ) =P((£( 7 r) ! £(a)) eC k ) 

= P (jJJnUl [nVTj g C fc } n {Aq n B n , fc » 

_ P (jjJnU] , [nVD e C k } n P n , fc ) + Q(P(P° |fc )) 

1 - o(P(p° ntk )) 

= P((U,V) GC k )+Q(k 2 /n) 

l-0(k 2 /n) 
= Q k + 0(k 2 /n), 

where Q k = P ((U, V) gCj). Let us write P n = P(tt < a). Using Q and the last estimate, we obtain then 

Pn > Q k P n -k (1 + 0{k 2 /n)) = Q k P n -k exp (0(k 2 /n)) , n > k. 
Iterating this inequality times gives 



Pi > Q k l/ki P^ln/klk exp 



[n/fej-l 
3=0 



n — jfc 
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Since the sum in the exponent is of order 0(k 2 logn), we get 



liminf yP n > \[Q k , Vfc > 1. 



Thus 



lim inf \/ P n > sup \J Q k 



k 



Therefore, for each k and e £ (0, y/Qk) , we have 



(5) 



Next 



Lemma 2 As a function of k, Q k is supermultiplicative, i.e. Qk t +k 2 ^ Qk x Qk 2 f or a ^ k\,k% > 1. Conse- 
quently there exists lim^oo y/Q k , and moreover 



Thus we expect that our lower bound would probably improve as k increases. 

PROOF. Q k is the probability of the event E k = {(U<*> , VW) £ C k }; here U<*> := {Ui}i<i< k , V« := 
{Vi}\<i<k- Explicitly, for each j < k and each c £ [0, 1], the number of l/i, . . . ,Uj not exceeding c is at 
most the number of V\,...,Vj not exceeding c. So Qk 1 +k 2 = P( E k 1 +k 2 )i Qfei = P{ E ki), while Q fc2 = 
P(Ek 2 ) = P(E% )■ Here the event means that for each j* < and each c G [0, 1], the number of U, 
i = ki + 1, . . . , kx + j*, not exceeding c is at most the number of Vi, i = k% + 1, . . . , k± + j*, not exceeding c. 
The events E kl and E^ 2 are independent. Consider the intersection of E kl and E% . There are two cases: 

1) j < fei- Then the number of U, i < j not exceeding c is at most the number of Vi, i < j not exceeding 
c, as E kl holds. 

2) k\ < j < k\ + ki- Then the number of Ui, i < j, not exceeding c is at most the number of Vi, i < k% 
not exceeding c (as E kl holds), plus the number of Vi, ki < i < j, not exceeding c (as E^ 2 holds). The 
total number of these Vi is the number of all Vi, i < j, that are at most c, c £ [0, 1]. 

So E kl +k 2 12 E kl n E^ 2 , and we get Q kl +k 2 > Qk x Qk 2 - The rest of the statement follows from a well-known 
result about super(sub)multiplicative sequences (see Polya and Szego [201 ) - 

Given 1 < j < i < k, let t/jj- (V^j resp.) denote the j-th element in the increasing rearrangement of 
U\, . . . ,Ui (Vi, . . . , Vi resp.). Then, to put it another way, Q k is the probability that the k Ehresmann condi- 
tions are met by the independent /c-dimensional random vectors U and V, both of which have independent 
entries. That is, we check Uj > V\j for each 1 < j < i < k by performing element-wise comparisons in the 
following tableaux: 




= sup 

k>l 




Uf.,1 £4,2 £4,3 



u kyk 



v kA v K2 v ky3 



v kyk 



U3.1 U3.2 C^3,3 

U2.1 U2.2 



^3,1 ^3,2 V3 ( 3 
V2.I V2,2 
^1,1 



What's left is to explain how we determined a = 0.70879. 
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It should be clear that whether or not (U( fe ) , V( fe )) is in C k depends only on the size ordering of 
Ui, . . . , Uk, Vi, . . . , Vfe. There are (2fc)! possible orderings, all being equally likely. Thus Qk = Nk/{2k)\. 
Since the best constant in the lower exponential bound is probably lim^oo \/Qk, our task was to compute 
Nk for k as large as our computer could handle. ("Probably", because we do not know for certain that \fQk~ 
increases with k.) 

Here is how N k was tabulated. Recursively, suppose we have determined all N k -i orderings of x\, . . . , 
Xk-i, 2/1, • ■ • , J/fe-i such that (x^ 0-1 ' , y^ -1 )) G Ck-i- Each such ordering can be assigned a 2(fc — l)-long 
sequence of O's and l's, O's for XiS and l's for t/j's, 1 < i, j < k — 1. Each such sequence meets the ballot- 
theorem condition: as we read it from left to right the number of l's never falls below the number of O's. 
We also record the multiplicity of each sequence, which is the number of times it is encountered in the list 
of all N k _i orderings. The knowledge of all 2(fc — l)-long ballot-sequences together with their multiplicities 
is all we need to compile the list of all 2fc-long ballot-sequences with their respective multiplicities. 

For k = 1, there is only one ballot-sequence to consider, namely 10, and its multiplicity is 1. So ^ = 1, 
and 

Qi = l/2L 

Passing to k = 2, we must count the number of ways to insert 1 and into 10 so that we get a 4- long ballot- 
sequence of two O's and two l's. Inserting 1 at the beginning, giving 110, we can insert into positions 2, 3 
or 4, producing three ballot-sequences 

1010, 1100, 1100, 

respectively. (Inserting into position 1 would have resulted in 0110 which is not a ballot-sequence.) 
Similarly, inserting 1 into position 2, we get 110, and inserting under the ballot condition gives three 
ballot-sequences 

1010, 1100, 1100. 

Finally, inserting 1 at the end, giving 101, we can only insert at the end, obtaining one ballot-sequence 

1010. 

Hence, starting from the ballot-sequence 10 of multiplicity 1, we have obtained two 4-long ballot-sequences, 
1010 of multiplicity 3 and 1100 of multiplicity 4. Therefore N 2 = 3 + 4 = 7, and 

Q 2 = 7/4!. 

Pass to k = 3. Sequentially we insert 1 in each of 5 positions in the ballot-sequence 1010, and then determine 
all positions for the new which would result in a 6-long ballot-sequence. While doing this we keep track of 
how many times each 6-long ballot-sequence is encountered. Multiplying these numbers by 3, the multiplicity 
of 1010, we obtain a list of 6-long ballot-sequences spawned by 1010 with the number of their occurrences. 
We do the same with the second sequence 1100. Adding the numbers of occurrences of each 6-long ballot- 
sequence for 1010 and 1100, we arrive at the following list of five 6-long ballot-sequences with their respective 
multiplicities: 

111000 : 36, 
110100 : 32, 
110010 : 24, 
101100 : 24, 
101010 : 19. 



Therefore N 3 = 36 + 32 + 24 + 24 + 19 = 135, and 
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Q 3 = 135/6!. 

The younger coauthor wrote a computer program for this algorithm. Pushed to its limit, the computer 
delivered the following table: 



k 


N k = (2k)\Q k 


Q k = N k /{2k)\ 


¥Qk 


1 


1 


0.50000 . . . 


0.50000 . . . 


2 


7 


0.29166... 


0.54006 . . . 


3 


135 


0.18750... 


0.57235... 


4 


5193 


0.12879... 


0.59906 . . . 


5 


336825 


0.09281 . . . 


0.62162... 


6 


33229775 


0.06937 . . . 


0.64101... 


7 


4651153871 


0.05335 . . . 


0.65790... 


8 


878527273745 


0.04198... 


0.67280... 


9 


215641280371953 


0.03368 . . . 


0.68608 . . . 


10 


66791817776602071 


0.02745 . . . 


0.69800 . . . 


11 


25497938851324213335 


0.02268... 


0.70879 . . . 



Combining (JSJ) and the value of S/Qn in this table, we see that for each e > 0, 

Pa = fi (( VOlT - e) n ) = ft ((0.708... - e) n ) . □ 

The numbers tyQ k increase steadily for k < 12, so at this moment we would not rule out the tantalizing 
possibility that \/Q k — > 1 as k — > oo. Determination of the actual limit is a challenging open problem. The 
proof just given only involves mention of the (0, l)-matrix criterion, but it was the Ehresmann criterion that 
actually inspired our initial insights. 

Weak Bruhat Order, preliminaries. 

Recall that tt precedes a in the weak Bruhat order (tt ^ a) if and only if there is a chain a = lo\ — > ■ ■ ■ — ► 
lu s = tt where each u>t is a simple reduction of ujt-i, i-e. obtained by transposing two adjacent elements 
u>t-i{i), Wt-i(i + 1) such that u>t-i(i) > u>t-i(i + !)• Clearly the weak Bruhat order is more restrictive, so 
that tt < a impies tt < a. In particular, P(ir < a) < P(n < a), hence (Theorem^ P(ir < a) = 0(n~ 2 ). 
We will show that, in fact, this probability is exponentially small. The proof is based on an inversion set 
criterion for <a implicit in ^ pp. 135-139]. 

Lemma 3 Given lu G S n , introduce the set of non-inversions of u: 

E(u):={(i,j) : i<j, u- 1 (i) < lj- 1 (j)} ■ 
i" ^ o~ if and only if E (tt) D E (a). 

Proof. Assume 7r ^< a . Then there exists a chain of simple reductions u>t, 1 < t < s, connecting a = lo\ 
and 7r = lo s . By the definition of a simple reduction, for each t > 1 there is i = it < n such that E(uj t ) = 
E(uj t -i) U {(uj t (i),u) t (i + 1))}, where u> t (i) = w t -i{i + 1), oj t (i + 1) = w t _i(i), and u t -i(i) > w t _i(i + 1). So 
the set E(u> t ) increases with t, hence E(n) 2 E(a). 

Conversely, suppose E (tt) D E (a). Since a permutation u) is uniquely determined by its E(u>), we may 
assume E (tt) D E (a). 
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Claim If E (tt) D E(a), then there exists u < v < n such that (v,u) is an adjacent inversion of cr, but 
(u,v) e E(tt). 

Assuming validity of the claim, we ascertain existence of an adjacent inversion (v, u) in a with (it, v) £ E (tt). 
Interchanging the adjacent elements u and v in a = lo\, we obtain a simple reduction LU2, with E(u>x) C 
E(u>2) C E(ir). If E(ix>2) — E(tt) then = 7T, and we stop. Otherwise we determine W3, a simple reduction of 
UJ2, with E(u>2) C E(lj^) C i?(7r) and so on. Eventually we determine a chain of simple reductions connecting 
a and 7r, which proves that tt < a . 

Proof of Claim. The claim is obvious for n = 1,2. Assume inductively that the claim holds for 
permutations of length n — 1 > 2. Let 7r, a £ S n and E(ir) D E(a). As in the proof of Theorem ^ let 
(-ni^) — 7T — 1 (tt.) , l n (cr) = er -1 (n), and tt* , a* are obtained by deletion of n from 7r and a. Since E (tt) D (cr), 
we have (-7T*) 3 i? (c*). Suppose first that (77*) = S (cr*). Then 7r* = cr*, and as E (n) D E (cr), we must 
have € n (7r) > £ n (a), i.e. ^ n (cr) < Setting v — n and u = cr (^„(cr) + 1), we obtain an adjacent inversion 
(i>, u) in cr with (u, v) £ E (tt). 

Alternatively, E (tt*) 2 E(tj*). By inductive hypothesis, there exists u < v < n — 1 such that (v,u) is 
an adjacent inversion of cr*, but (u, v) £ E(jr*). Now insert n back into 7r*,cr*, recovering n and cr. If n 
sits to the right of u or to the left of v in cr, then (w,u) is still an adjacent inversion of cr. Otherwise n is 
sandwiched between v on the left and u on the right. Therefore (n, u) is an adjacent inversion in a. On the 
other hand (v, n) £ E (a), so since E (tt) D E (cr), we have (w, n) £ E (tt) also. Hence, the triple (u, v, n) are 
in exactly this order (not necessarily adjacent) in tt. Therefore the adjacent inversion (n, u) in a is such that 
(it, n) £ E (tt), and this proves the inductive step. D 

Denote by Q the permutation lu reversed in rank. For example, with tu = 13254 we have u) = 53412. Then 
it is easy to see that 



E (tt) D E (cr) 
By Lemma these statements are equivalent to 



E (tt ) C E (a) . 




We immediately obtain the following corollary to Lemma [21 



Corollary 1 For u £ S n , define 



Ei H := {j < i : (j,i) £ E (u>)} , 1 < i < n. 



Then 



n 



and consequently 




□ 



Next. 
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Lemma 4 For tt,o~ G S n selected independently and uniformly at random, let P* :— P (tt -< cr). P* is 
submultiplicative, i.e. for all n\,n 2 > 1 

P* <C P* P* 
ni+n2 — Hi n-2 ' 

Consequently there exists liirin^oo \/P* = inf„>i y/P£. 

Proof. Let tt, cr be two permutations of [ni + na]- Then tt ^ cr if and only if 

#i (tt) 2Ei(a), 1 < i < ni + n 2 . 
Using these conditions for i < n\, we see that 

7r[l,2,...,n 1 ]^ff[l,2,... > ni]. 

Here tt [1,2, . . . , m], say, is what is left of the permutation tt when the elements ni + 1, . . . , n\ + n 2 are 
deleted. 

Likewise, 7r X cr if and only if 

^ (tt) CEi(a), 1 < i < n x + n 2 . 
Using these conditions for i < n 2 , we see that 

7r [n x + 1, . . . , rii + n 2 ] -< a [rii + 1, . . . ,ni + n 2 ] ■ 
Now, since tt and cr are uniformly random and mutually independent, so are the four permutations 

7r[l,...,m], Tr[n 1 + l,...,n 1 +n 2 ], a [1, . . . , m] , cr [ni + 1, . . . , rii + n 2 ] . 

Hence, 

P (tt r< cr) < P (tt [1, . . . , ni] d cr [1, . . . , m]) ■ P (tt [nx + 1, . . . , m + n 2 ] r< <r [m + 1, . . . , m + n 2 ]) , 
so that 

p* <^ p* p* n 
Proof of Theorem El upper bound. 

We will show that for each e > 0, 

P* = 0((<3 + eT), 

where 



The proof of this upper bound for P* parallels the proof of the lower bound for P„ in Theorem ^ As in 
that proof, given k > 1, let ir k * and cr fc * be obtained by deletion of the elements n, . . . , n — k + 1 from tt and 
cr, and let £i(ir) = 7r _1 (i), ii(cr) = cr _1 (i), n — k + 1 < i < n. In the notations of the proof of Lemma^J 
iT k * = 7r[l, . . . , n — k] and a k * = cr[l, . . . , n — k], and we saw that TT k * ^ a k * if tt -< a. Our task is to find the 
conditions these £,(-)'s must satisfy if tt -< a holds. 
To start, notice that 
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7T X a => |£ n (7r)| > |£„(a)| <=> 4(0 > 4(<r)- 

Next 

7T X a => 7T* X a* => 4_!(7r) > 4-1 (O > 4-1 (O > 4-i(<t) - 1, 

as deletion of n from 7r, ct decreases the location of n — 1 in each permutation by at most one. In general, 
for < j < k we get 

7T < a tt 3 ** d > - J. 

So, introducing £(tt) = {4-i+i(n")}i<i<* and £(cr) = {4-i+i(c)}i<i<fc, 

d ^} C {(£(7r),/(tr) G (6) 
5 fc := {(x,y) G M 2fe : x, > Uj - j + 1, 1 < j < k} . 

In addition, on {-7T ^< er} every pair of elements, which forms an inversion in ir, also forms an inversion in a. 
Applying this to the elements n — k + 1, . . . , n, we have then 

{w^a}C{(£(n),£(a))€T k }, (7) 
% := {(x, y) G R 2k : VI < i < j < k , x, < x, } ■, Vi < Vj }. 

Combining © and 10, we get 

{tt < a} C {(£(ir),£(a)) £ S k n %} n {ir k * * a k *} . 
So, since the two events on the right are independent, 

P'n < P((tWA°)) %)P* n _ k . (8) 

It remains to estimate P((£(ir), £(a)) G Sfc D Tfc). As in the proof of Theorem^ (lower bound), we observe 
that (£(tt) , £(&)) has the same distribution as (|~7iU~|, |~nV~|), conditioned on 

An, k n B n>k = { \nU{\ \nU k ] } n { \nV{\ \nV k ]}. 

Here U\, . . . , U k , V\, . . . , V k are independent [0, l]-uniforms. Then 

P((fM HtrWc: <? nTl P({(WJnVl)e5 fc nr fc }nc n , fc ) 

e n ifcj = — — — r , L nt k = Ai,fc n 

Introduce the event 2? n! fc on which 

minjmin |Z7j — min |V£ — Vj\, min|Z7j — Vj\, fc _1 min \Uj — Vj\} > 1/n. 
Certainly 2?„.fc C C„ fc and, thanks to the factor 1/fc by min^ \Uj — Vj\, on T> n ^ k 

\nUj] > \nVj] - j + l^Uj>Vj - k/n=> Uj > V r 

Therefore, on T> ntk , 

( [nUl , \nV] )es k n%=^ (U, V) G s k n r fe , 
5 fc := {(x,y) G R 2k : x d > Vj , 1 < j < k}. 
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Clearly S k fl Tk is a cone-shaped subset of R 2fe . In addition, P(T>^ k ) = 0(k 2 /n). Hence 

PI £(7r ),l (a)) G Ofc n Tfe) < = 

l-0(P(D° ik )) 

= Q* k (l + 0(k 2 /n)), Ql :=P((U,V) G5 fc nT fc ). 

This and JHJ) imply 

P:<QtP^ k exp(0(k 2 /nj). 
Hence, as in the proof of Theorem ^ (lower bound), 

limsup tfP* < {/Q*, Vfc > 1, 

and so 

Pn=O(($/0l + e) n ), k>l,e>0. (9) 
Furthermore, from the definition of Q k , it follows directly that Q* k is submultiplicative, i.e. 

Ql 1+k2 < Q* kl Ql 2 , k u k 2 >i. 

Therefore ([20] again) 

lim VO* = inf {/Ql . 

k — >oc k>l 

So the further we can push tabulation of Q* k , the better our exponential upper bound for P* would probably 
be. ("Probably" , because we do not have a proof that \fQ% decreases with k.) 

As in the case of Q k , Q k = N k /{2k)\. Here, by the definition of the sets S k and T k , N k is the total 
number of ways to order x\, . . . , x k , y\, . . . , y k so that two conditions are met: (1) for each j, Xj is to the 
right of yj] (2) for all i < j, if Xi is to the left of Xj then y\ is to the left of yj. 

It is instructive first to evaluate N k by hand for k = 1,2. N± — 1 as there is only one sequence, yiX\, 
meeting the conditions (1), (2). Passing to iV|, we must decide how to insert yi and x 2 into the sequence 
y\X\ in compliance with conditions (1), (2). First of all, y 2 has to precede x 2 - If we insert x 2 at the beginning 
of yxXi, giving x 2 y±xi, then we can only insert y 2 at the beginning of this triple, giving 



y 2 x 2 yxxi. 

Alternatively, inserting x 2 in the middle of y\X\, we have 2 possibilities for insertion of y 2 , and we get two 
admissible orderings, 



y 2 y\x 2 x\, yiy 2 x 2 xi. 

Finally, insertion of x 2 at the end of y\X\ brings the condition (2) into play as we now have x\ preceding x 2 , 
and so j/i must precede y 2 . Consequently, we get two admissible orderings, 



yiy 2 xix 2l yixiy 2 x 2 . 

Hence N 2 =1 + 2 + 2 = 5. Easy so far! However, passing to k = 3 is considerably more time-consuming than 
it was for computation of N3 in the proof of the lower bound in Theorem^ There, once we had determined 
the N 2 admissible orderings, we could afford not to keep track of relative orderings of X\, . . . , x k -\, and of 
yi, ... , ?/fe-i, whence the coding by l's and O's. All we needed for passing from k — 1 to k was the list of all 
binary ballot-sequences of length 2(fc — 1) together with their multiplicities. Here the nature of the conditions 
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(1), (2) does not allow lumping various sequences together, and we have to preserve the information of relative 
orderings of x's, and relative orderings of j/'s. This substantial complication seriously inhibits the computer's 
ability to compute JV£ for k as large as in the case of Nk ■ 

To get a feeling for how sharply the amount of computation increases for k = 3, let us consider one of 
the N£ = 5 admissible sequences, namely 2/2X22/1X1. As above, we write down all possible ways to insert 2/3 
and X3 into this sequence so that (1) and (2) hold. Doing this, we produce the 10 sequences: 



2/3X32/2X22/1X1, 
^22/3X3X22/1X1, 
2/2X22/3X32/1X1, 

2/2^22/32/1^3^1, 

2/2X22/12/3x1x3, 



2/32/2^3X22/1X1, 
2/22/3X2X32/1X1, 
2/22/3X22/1X3X1, 
2/2X22/12/3X3X1, 
2/2X22/1X12/3X3. 



We treat similarly the other four sequences from the k — 2 case, eventually arriving at N£ = 55. We wouldn't 
even think of computing N£ by hand. 

Once again the computer programming to the rescue! Here is the table produced by the computer after 
a substantial running time: 



k 


N* = (2k)\Ql 


Ql = N*/(2k)\ 




1 


1 


0.50000... 


0.50000 . . . 


2 


5 


0.20833... 


0.45643 . . . 


3 


55 


0.07638 . . . 


0.42430 . . . 


4 


1023 


0.02537... 


0.39910... 


5 


28207 


0.00777... 


0.37854 . . . 


6 


1065317 


0.00222 . . . 


0.36129... 



Using © for k — 6, we get for each e > 

p: = ((Voi+^) n ) = ((o-36i...+en. □ 

Proof of Theorem [21 lower bound. 

To bound P (tt ■< a) from below we will use the criterion 

7T < a Ei (tt) D E l (a) , Vi < n. 

First of all, 

Lemma 5 Let i G [n], B C [i — 1] ( [0] = $). Ifir£S n is chosen uniformly at random, then 

P (Ei (tt) DB)= ' 



B\ + 1 

PROOF. By the definition of Ei(ir), 

{EM DB} = {jr-^j) < n^d), Vj G B}. 
It remains to observe that n^ 1 is also uniformly random. D 

Lemma |5] implies the following key statement: 
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Lemma 6 Let it, a G S n be selected independently and uniformly at random. Then, for i G [n], 



PROOF. By Lemma 



1 

P(E i (w)DE i (a))=H(i)/i, ff(i):=V-. 

7 



F(£,W2£,M)= J! P(E i (7r)2B)P(E i (a) = B) 

BC[i-l] 

P(E i (a) = B) 



E 



LB +1 



T7 



1 _#(*) D 
f-r( i fi + H « 



Note. In the second to last equality, we have used the fact that |-Ej(<r)| is distributed uniformly on 
{0, 1, . . . , i — 1}. In addition, |23i(<r)|,. . . , |-E„(ct)| are independent, a property we will use later. For com- 
pleteness, here is a bijective proof of these facts. By induction, the numbers |J3j(<r)|, i < t, determine 
uniquely the relative ordering of elements 1, . . . , t in the permutation a. Hence the numbers |i?i(cr)|, i G [n], 
determine cr uniquely. Since the range of | £7j (cr) | is the set {0, . — 1} of cardinality i, and \S n \ = n\, it 
follows that the numbers |-Ej(er)|, i£ [n], are uniformly distributed, and independent of each other. 

Needless to say we are interested in P(tt :< cr) = P (n" =1 {£'i(7r) D Ei(cr)}). Fortunately, the events 
{Ei(jr) D Ei(a)} turn out to be positively correlated, and the product of the marginals P(Ei(n) 2 Ei(a)) 
bounds that probability from below. 

Theorem 3 Let tt,ct G S n be selected independently and uniformly at random. Then 

n n jj , 

p(tt d a) > nm w 2 E i w) = n ~i- 

i=l i=l 
Proof. First notice that, conditioning on cr and using the independence of it and cr, 

P{E l {n) 2Ei{o-),Vi<n) 

= E[P{E i {K)DE i (a),Vi<n\a)}=E[P(E i (- ! r)DB i ,Vi<n)\ Bi = Ei(a) ] . 

So our task is to bound P (Ei (n) D Bi, \fi < n), where these Si's inherit the following property from the 
^ (a)'s: 

i G Ej (cr) and j G E k (cr) =>• i G E k (cr) . 
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Lemma 7 Let n > 1 be an integer, and let Bi C [n], i = 1, . . . ,n, &e smc/i t/iai 

i £ Bi and i G Bj, j G £?fc i G £>fc, V z, j, fc G [n]. 
Then, for ir € S n selected uniformly at random, 



P (Ei (tt) D S ia V* < n) > f[ T^TTY- 



Proof of Lemma El Notice upfront that UiBi ^ [n]. Otherwise there would exist ii,.. .,i s such that 
it G Sj {+1 , 1 < t < s, (i s +i — and - using repeatedly the property of the sets Bi - we would get that, 
say, ii G Bi 2 and 12 G S^, hence «2 G S^ 2 ; contradiction. 

Let U±, . . . , U n be independent uniform- [0, 1] random variables. Let a random permutation u be defined 

by 

w{i) = k •<=>■ /7i is fc th smallest amongst Ui, . . . ,U n . 
Clearly u> is distributed uniformly, and then so is tt := tu^ 1 . With tt so defined, we obtain 

{E, (tt) D Bi, Vi < n} = {tt- 1 (i) > tt' 1 (j) , V j eB z ,i<n} = {U > U jt Vj G B t , i < n} . 
Hence, the probability in question equals 

P{U i >U ji Vj€B i ,i<n). 
We write this probability as the ra-dimensional integral 

P (Ui > £/j, V j G Bi, i<n) = J ■ ■ ■ J dxi ■ ■ ■ dx n , 

D 

D = {(xi, . . . ,x n ) e [0, l] n : Xi > Xj, V j G Bi, i < n} . 

Since UiBi =/= [n], we can choose an index k G [n] such that k ^ Bi for all i. Then we may rewrite the integral 
above as 



Jo [ J ■ ■ ■ J dx! ■ ■ ■ dx k -idx k +i ■ ■ ■ dx n dx k , 

\D(x k ) J 

(xk) = . . . ,Xk-i,Xk+i, ■ ■ • , x n ) G [0, : x t > x 3 , V j G Bi, i < raj . 



D 



On D (xk), the only inequalities involving Xk are of the form Xk > Xj, j G B k . This suggests scaling those Xj 
by Xk, i.e. introducing new variables tj := Xj/xk, so that tj G [0, 1], j G Bk- To keep notation uniform, let 
us also replace the remaining Xi, i ^ Bk U {k}, with ti. Let 2) (xk) denote the integration region for the new 
variables U, i ^ k. Explicitly, the constraints Xj < Xk, j G Bk, become tj < 1, j G Bk- Obviously each listed 
constraint x a < Xb (a, b G Bk) is replaced, upon scaling, with t a < tb. We only rename the other variables, 
so every constraint x a < Xb (a, b (£ Bk) similarly becomes t a < tb. By the property of the sets Bi, there are 
no inequalities x a > Xb, a G Bk, b ^ Bk (since the presence of this inequality implies b G B a ). The only 
remaining inequalities are all of the type x a < Xb, a G Bk, b ^ Bk- In the new variables, such a constraint 
becomes Xkt a < tb, and it is certainly satified if t a < tb, as Xk < 1. Hence, !D (xf.) 2 D* , where 

D* := {(h, . . . ,t k ^,t k+l , . . . ,t n ) G [0,!]™- 1 : U > tj, Vj G B u i ? k} , 
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and D* does not depend on x^- Observing that the constraints that determine D* are those for D with the 
constraints < Xk, i S -Bfe, removed, we conclude that the innermost integral over D(xk) is bounded below 
by x k * P{Ui > Uj, V j G B il i ^ k). (x k k is the Jacobian of the linear transformation {2^}^ — > {ti}i^k-) 
Integrating with respect to Xk, we arrive at 



P(U, > U jt Vj € £i, i < n) > 



1 



P(Ui>Uj,Vj €Bi,i?k). 



(10) 



By induction on the number of sets Bi, with Lemma [S] providing basis of induction and I|1U|I - the inductive 
step, we get 



" 1 

P (Ui > U h Vj €B z ,i< n) > J[ j^ppy 
The rest is short. First, by Lemma [7| 

P(£7i(7r) D BiCtr), Vz < n) = £ [p(^(tt) 2 B,, Vz < 



□ 



l_Bi=_Bi(cr) 



> E 



n 



11 | J B i (a)| + l_ 

Since the cardinalities \E^ (a) | are independent, the last expected value equals 



i=l 



\Ei (a) I + 1 



IIItEtttI-II 



Hi). 



, 1 * — ' 7 + 



for the second to last equality see the proof of Lemma 



□ 



Note. Let V be a poset on [n], and put := {.?' e : j < i in P}. Bi U {i} is called the order ideal at i. 
By the properties of V, the Bi's satisfy the hypotheses of LemmaQ so letting e (V) denote the number of 
linear extensions of V we get 



P(Ei(7r)2B h V«<n) = 



I : w (i) > oj (j) , Vj G Bi, i <n}\ 



e(V) 



> 



- 1 



Thus we have proved 



Corollary 2 For a poset V with n elements, 

n 

e{V) > nl/jj <*(*)> :=|{i6P:j<imP}|. □ 



In a very special case of V , whose Hasse diagram is a forest of rooted trees with edges directed away from 
the roots, this simple bound is actually the value of e(V) (see JHI Section 5.1.4, ex. 20]). There exist bounds 
for the number of linear extensions in the case of the Boolean lattice (see Brightwell and Tetali 6 , Kleitman 
and Sha ^lj), but the lower bound in Corollary [2 seems to be new. 
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Numerics 

From computer-generated data we have collected, it appears that our O (n~ 2 ) upper bound correctly 
predicts the qualitative behavior of P (n < a). The data suggests that P (n < a) is of exact order n~' 2+(5 ^ for 
some 5 € [0.5, 1], which begs the question of how to improve on our current bound. Writing P n = P (ir < a) , 
below is a graph (based on this numerical experimentation) exhibiting convergence to the exponent —a in 
the asymptotic equation P n <~ cn~ a , c > a constant, and —a appears to be near —2.5: 



10 



Value of « 

30 50 70 90 110 



-1 



-1.5 



5 

S 



-2 



-2.5 



Here is a portion of the accompanying table used to generate this graph: 



n 


Rn 


Estimate of P n s= 


_ Rn 


Estimate of ln(P„)/ In n 


10 


61589126 


0.0615891... 




-1.21049... 


30 


1892634 


0.0018926 . . . 




-1.84340... 


50 


233915 


0.0002339 . . . 




-2.13714... 


70 


50468 


0.0000504 . . . 




-2.32886... 


90 


14686 


0.0000146 . . . 




-2.47313... 


110 


5174 


0.0000051 . . . 




-2.58949... 



Here, R n is the number of pairs (n, a) out of 10 9 randomly-generated pairs such that we had ir < a. We 
have also utilized the computer to find the actual probability P n for n = 1, 2, . . . , 9. Below is a table of these 
true proportions: 



n 


(n\) 2 P n 


Pn 


1 


1 


1.00000... 


2 


3 


0.75000 . . . 


3 


19 


0.52777... 


4 


213 


0.36979 . . . 


5 


3781 


0.26256... 


6 


98407 


0.18982... 


7 


3550919 


0.13979... 


8 


170288585 


0.10474... 


9 


10501351657 


0.07974 . . . 
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Concerning the weak Bruhat order, computer-generated data suggests that P (tt < a) is of exact order 
(0.3)™. So our current upper bound O((0.362) n ) is a qualitative match for P {it -< a), but it appears that 
improvements are possible here also. Writing P* = P (tt ^ a), below is a graph (based on our numerical 
experiments) exhibiting convergence to the ratio p in the asymptotic equation P* <~ cp n , c > a constant, 
and p appears to be near 0.3: 



l.i 



0.9 



S • t 0.7 

U Oh 

.a 



Z 0.5 



0.3 



Value of n 
1 



10 



13 



Here is a portion of the accompanying table used to generate this graph: 



n 


H n 


Estimate of P* ? 


" T7F 


Estimate of P*/P*_i 


10 


1538639 


0.0015386... 




0.368718... 


11 


541488 


0.0005414 . . . 




0.351926... 


12 


184273 


0.0001842 . . . 




0.340308 . . . 


13 


59917 


0.0000599 . . . 




0.325153... 


14 


18721 


0.0000187... 




0.312448 . . . 


15 


5714 


0.0000057... 




0.305218... 


16 


1724 


0.0000017... 




0.301715... 



Here, R* n is defined analogously to R n above. Below is a table with the true proportion P* for n = 1, 2, . . . , 9: 



n 


(ni) 2 p r ; 


P* 


1 


1 


1.00000. . . 


2 


3 


0.75000. . . 


3 


17 


0.47222.. . 


4 


151 


0.26215.. . 


5 


1899 


0.13187.. . 


6 


31711 


0.06117.. . 


7 


672697 


0.02648. . . 


8 


17551323 


0.01079.. . 


9 


549500451 


0.00417.. . 
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Surprisingly, our Theorem [3 lower bound for P* is quite good for these smallish values of n: 



n 




l 

2 
3 
4 
5 
6 
7 
8 
9 



1569.8. 
23075.9. 
418828.3. 
9106523.1. 
231858583.9. 



1.0. 
3.0. 
16.5. 
137.5 . 



1.00000. 
0.75000 . 
0.45833 . 
0.23871 . 
0.10901 . 
0.04451 . 
0.01648 . 
0.00560 . 
0.00176. 
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